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We develop a rigorous mathematically consistent description of PT symmetric optical systems 
by using second quantization. We demonstrate the possibility of significant spontaneous generation 
of photons in PT symmetric systems. Further we show the emergence of Hanbury-Brown Twiss 
(BIBT) correlations in spontaneous generation. We show that the spontaneous generation determines 
decisively the nonclassical nature of fields in PT symmetric systems. Our work can be applied to 
other systems like plasmonic structure where losses are compensated by gain mechanisms. 
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The PT symmetric Hamiltonians have been extensively 
studied in quantum mechanics [iHl]. Generally the PT 
symmetric Hamiltonians involve complex potentials V{x) 
with the property that V{—x) — V*{x). It has been well 
known in quantum mechanics that a complex potential 
implies a source or sink of energy. The well known conti- 
nuity equation has a nonconservative term which is pro- 
portional to V{x) — V*{x) = V{x) — V{—x). Therefore a 
consistent theoretical description should account for such 
sources or sinks of energy. In optical systems the refrac- 
tive index is like a potential Q and since one can engi- 
neer different refractive indices, one effectively generates 
a very wide class of potentials leading to optical real- 
ization of PT symmetric potentials as demonstrated by 
Ruter et al0-@|- In optical realization of PT symmetric 
Hamiltonians one considers overall potential to be com- 
posed of two layers — layer A lying in the domain say 
— 1 < a; < and a layer B in the domain in < a; < +1. 
The two layers are characterized with complex optical 
susceptibilities x related to complex refractive index via 
n'^ — 1 + 4:Trx which are homogeneous and which have 
properties Re[xA] = R.e[xB] and Imfx^] = -Imtxs]- 
The layer A could be a gain medium and then the layer B 
would be an absorbing medium. Since the medium A is 
a gain medium it must be pumped by an external source. 
The medium B dissipates energy to an external environ- 
ment. We must treat consistently these dissipative or 
gain environments, and this shows how PT symmetric 
optical systems can lead to significant spontaneous gen- 
eration of photons. The possibility of spontaneous gen- 
eration has been missed so far and would be important 
when examining the transmission of single photons and 
entangled photons through PT symmetric structures [lol- 
12|. The transmission of single photons and in general 
quantum states of light through optical systems with 
Im[x^] — Im[xs] = is currently being studied exten- 
sively, particularly in quantum information science and 
in problems like quantum random walks 14 1 . 

Consider the PT symmetric optical system shown in 
Fig. [H The system consists of two single mode waveg- 
uides which are evanescently coupled. The waveguide A 
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Input states— vacuum, single 
photon or entangled states 

FIG. 1: The PT symmetric optical system consisting of 
two waveguides coupled evanescently; the waveguide 
A(B) is filled with gain(lossy) medium. 



is filled with a gain medium whereas the waveguide B 
is filled with a lossy medium. The PT symmetry would 
require that the gain in A must be equal to absorption 
in B. Let a and /? be the classical field amplitudes in 
the two waveguides. These amplitudes satisfy the basic 
equations 
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(1) 



We would now discuss the equations for the quantized 
fields in A and B. Let a, b {a\ h^) be the annihilation 
(creation) operators for the fields in A and B. In order to 
obtain the quantized equations, we note that any gain or 
lossy medium leads to quantum fluctuations in the fields. 
The basic equations can be derived from considerations 
of the microscopic Hamiltonian for the radiation fields, 
pumping and loss mechanisms. Calculations show that 
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Ga - iJb + fa, 



-Gb - Ua + fb 
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Here fa and /& are the quantum Langevin forces with 
zero mean and with the properties 

(/t(z)/„(z')> = 2G6{z - z'), {fa{z)fl{z')) = 0, 
{fl{z)h{z')) = 0, {h{z)fl{z')) = 2G5{z - z'), 

(3) 

Further the fluctuating forces are Gaussian. Note the dif- 
ference in the nature of Langevin forces in the equations 
for the amphfying and absorbing medium. These quan- 
tum mechanical fluctuating forces in general would con- 
tribute to the output of the waveguides and significantly 
modify the outputs. The quantum Langevin equations 
([2]) can be integrated over the length of the waveguide. 
We write the solution as 

where 

i = Jz^ ,^§. (5) 

The matrix M has the eigenvalues ±i-\/l — if g = 
G/J < 1, ±-\/ — 1 if g > 1. The intensities at the out- 
put I = Jz are obtained by calculating la = {a^{l)a{l)), 
h = (&^(0^(0) ^ variety of input states. The in- 
put states could be single photon states like |1,0), |0, 1) 
and even a NOON state like (|0,2) -I- |2,0))/v^. We 
next give explicit results for a variety of input quantum 
states in the PT symmetric optical system. Note that the 
biorthogonal modes of the PT symmetric optical system 
are given by the eigenstates of M. In view of the presence 
of quantum Langevin terms in Eq.Q the eigenmodes of 
the system are strongly correlated. 



Spontaneous Generation of Photons in PT 
Symmetric Optical Systems 

Let us first consider the most intriguing possibility of 
spontaneous generation in PT symmetric systems. In 
this case the input state of photons is |0, 0). Further the 
quantum noise in the absorber does not contribute to 
the normally ordered moments and thus only fa terms 
contribute to spontaneous generation. Using Eq.(j31) and 
(|3]) we can show that the spontaneous generation is given 
by 

Ia{l) ^ Sail) ^ 2g [ dl'\Kaain\^ 

% (6) 

h{l) = Sb{l) = 2g / Al'\K^a{l')\\ 
Jo 

where Kij = (e*^')y . We show in Fig. [51 the nature 
of the spontaneous generation for different values of the 



gain both below and above threshold g = I. The thresh- 
old is defined by the condition that the eigenvalues of the 
matrix M change from pure imaginary to pure real val- 
ues. We note from the Fig. [5] that the generation in the 
waveguide B starts later as initially the generation has 
to take place in the waveguide A and this must tunnel 
into the waveguide B to see the generation in the waveg- 
uide B. For large g the generation in the waveguide A 
dominates over that in B 




(a) (b) 



FIG. 2: Spontaneous generation of the radiation in PT 
symmetric optical systems. The solid and dashed curves 
give respectively the output of the waveguides A and B. 
Notice the dashed curves are multiplied by a factor of 

10. 



Vacuum Induced Hanbury-Brown Twissfl5l-[T7t 
Cross Correlation 

We next examine the cross correlation between the out- 
puts of the two waveguides. We expect such a cross cor- 
relation to be nonzero as the generation in B is due to 
generation in A and due to tunneling effects. Let us de- 
fine a correlation coefficient as 

qil) = g<^'^~l, 
(2) ^ {aHl)bHl)ami)) (8) 
^ (at(0a(0>(6t(06(0>- 

The term {a'^ {l)b'^ {l)a{l)b{l)) is proportional to the prob- 
ability of simultaneously detecting one photon in the out- 
put of both A and B waveguides. This can be evaluated 
by using the Gaussian nature of the Langevin forces fa 
and /h which implies that the Heisenberg operators are 
also Gaussian with the provision that we need to keep 
track of the operator characteristics. It can be shown 
using the Gaussian property that 

{a^ (Ip {l)a{l)b{l)) 

= (at(0a(/))(6t(05(/)) + {a\l)b{l)) {b\l)a{l)) (9) 
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and hence 



9(0 = 
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The term {a^{l)b{l)) is obtained from Q as 
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The nature of the vacuum induced HBT cross correlation 
between the outputs of the two waveguides is shown in 
Fig. [21 The behavior of the HBT correlation is quite 
different for g < 1 and g > I. For g > 1, the correlation 
saturates to unity, whereas for 5 < 1, it tends to remain, 
in most cases, much smaller than unity. Thus for large 
gain, the correlations or the fluctuations are maximally 
Gaussian. 




FIG. 3: Vacuum induced Hanbury-Brown Twiss (HBT) 
cross correlation between the outputs of the two 
waveguides. 



Stimulated vs Spontaneous Generation in PT 
Symmetric Optical Systems 



Hence asymptotically negligibly small output is there 
from the waveguide B. In Fig. |4] we present explicit re- 
sults for the stimulated generation and the total genera- 
tion for the input single photon either in the waveguide 
A or in B. The stimulated part is same as from Eq.([T]) 
and hence labeled as classical in Fig. 2] The behavior is 
quite different for g < 1 and for g > 1. For 5 < 1 , we find 
that total generation of photons is dominated by spon- 
taneous generation. The situation is somewhat different 
for g > 1. However, spontaneous generation still remains 
quite important. 
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FIG. 4: The mean number of photons in the output of 

waveguides A (solid curves) and B (dashed curves) 
when the input state is |1,0) (panels a, b), |0, 1) (panels 
c, d). Panels a and c are for g = 0.9; panels b and d for 
5 = 1-5. 



We would now consider the stimulated generation by 
single photon input states and compare the stimulated 
generation with spontaneous generation. The stimulated 



generation comes from the terms e 



Ml 



/a(0) 



Eq.dH). 



For the input |1, 0) the stimulated generation is given by 



(12) 



whereas for the input single photon in the waveguide B, 
the stimulated generation is given by 



/f = |i^,,(OP, = |iffc,(/)|^ 



(13) 



Note that for large propagation distance and g > 1, the 
stimulated generation is determined by 
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\Kab{l)\^ ^ \Kbail)\' 
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(14) 



4g2 



Hanbury-Brown Twiss Correlations for Input Fields 
in a NOON State 



Lastly we examine the behavior of optical fields in PT 
symmetric systems when the input fields are in a NOON 
stated!] (|0,2) + |2,0))/%/2 which is an entangled state. 
For this state g^^^ has the least value zero. The mean 
intensity of emission has a behavior similar to that in 
case of the state |1, 0) 



(a^a) = \Kaa\'' 



{b^b) 



bb\ 



\Kab\'^ + Sa, 



Sb- 



(15) 



The probability of finding one photon in each waveguide 
is the sum of three different types of contributions: (i) 
stimulated ones \KaaKba + KabKbbl'^', (ii) purely from 
spontaneous emission SaSb + SabSba', (hi) cross terms in- 
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volving both spontaneous and stimulated emission 

Sa{\Kba\'' + \Kbb?) + Sb{\Kaa? + \Kab\'') 

+ {SabKl^Kaa + SabK^Kab + C.C.}. (16) 

We show in Fig. [5] the correlation g^^^ in case if we com- 
pletely ignore spontaneous generation, i.e. the result with 
Sa = Sb = Sab = 0, and compare with the full result — 
which includes all the three contributions as mentioned 
above. Fig. [5t and Fig. [5ji show that g^^-* is mostly 
less than unity if the spontaneous generation is ignored 
whereas the inclusion of spontaneous emission changes 
the character of g''^) completely (Fig. [5^ and Fig. [SId). 
Thus noninclusion of spontaneous emission would lead 
to erroneous conclusion that the fields in FT symmet- 
ric system are nonclassical whereas the correct behavior 
leads to nonclassicality only over a small range of propa- 
gation distance. 
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FIG. 5: HBT cross correlation between the two 
waveguides when the input is a NOON state 
(|2,0) + \0,2))/y/2. In panels a and b, we consider both 
spontaneous and stimulated generation; while in panels 
c and d, we consider only the stimulated one. 

In summary we have shown how important is the spon- 
taneous generation of photons in FT symmetric optical 
systems. We show significant Hanbury-Brown Twiss cor- 
relations in the generated photon fields. The sponta- 
neously generated photons strongly correlate eigenmodes 
and determine the quantum nature of the fields contain- 
ing single and entangled photons in such systems. There 
are many other important systems involving for exam- 
ple the lossless propagation in plasmonic structures [l9j| 



where losses are compensated by gain media. While deal- 
ing with quantum fields in such structures one has to 
supplement Maxwell equations with quantum Langevin 
terms of the type discussed in the current work. These 
ideas would also be applicable, say, in the context of mat- 
ter waves where a whole class of potentials including com- 
plex ones can be constructed by using off resonant laser 
fields^]. 
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